Charge and spin response functions for the Tomonaga model with quadratic 

dispersion and different interactions 
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We derive expressions for the charge and spin response function for the Tomonaga model with 
quadratic dispersion and arbitrary (but finite for zero momentum) interaction. For constant inter- 
action these expressions are analytic and for other types of interaction only a simple matrix has to 
be diagonalised. We use a truncated expansion in particle-hole states with and without inclusion 
of correlations in the ground state yielding an exact result for pure intra-band interaction. We also 
discuss the possibility of power-laws in the dynamic structure factor for the spinful and spinless 
model. 



I. INTRODUCTION 

Low dimensional systems of interacting electrons are 
an active field of research, with the interaction playing 
a more dominant role than in higher dimensions. In his 
famous 1950 work Tomonaga 1 found the exact solution 
of such a model system for one spatial dimension based 
on the following assumption: If the range of the interac- 
tion is much larger than the inter-particle distance, the 
quadratic energy dispersion can be linearised around the 
Fermi points. Tomonaga found collective excitations of 
the electronic density, i.e. plasmons as the low-lying ex- 
citations of such a system. The systematic extension of 
this model is nowadays called Tomonaga-Luttinger (TL) 
model and discussed in many reviews^ and books^ 

Quite recently, several authors attempted to go be- 
yond this approximation and to include effects of the 
non-linearity in the energy dispersion.— Of main 
interest was the shape of the dynamic structure factor 
(DSF) where the plasmons of the TL model show up as 
simple delta-peaks and changes are expected as broaden- 
ing of the peaki (corresponding to plasmon damping) or 
as the appearance of power-laws.^^ 
A possible starting point for a discussion could be given 
by the random phase approximation (RPA) which is 
nothing else but linearised time-dependent Hartree the- 
ory with the Hartree-Fock energies replaced by the non- 
interacting ones. This approximation becomes exact in 
the limit of linear dispersion, yielding the exact solution 
of the TL model. But it has been shown&& only recently 
as an insufficient starting point for a discussion of the 
above mentioned questions. An analytic solution for the 
special case of a constant interaction has been derived 
within a better approximation - the random phase ap- 
proximation with exchange (RPAE) corresponding to lin- 
earised time-dependent Hartree-Fock theory (TDHF). 

The linearization of the dispersion leads to two linear 
branches instead of a single quadratic one. Electrons on 
the right branch are usually called right-movers and those 
on the left branch left-movers. Accordingly, one can dis- 



tinguish interactions between electrons on the same and 
on different branches. In this work we make explicit use 
of this distinction between inter-band and intra-band in- 
teractions. The shape of the interaction under consider- 
ation is arbitrary with only one constraint necessary for 
Luttinger liquid phenomenology: the Fourier components 
of the interaction are finite and only non-zero for small 
momenta. 3 After introducing response functions in gen- 
eral in Sec. [ill we focus on pure intra-band interactions 
in Sec. IIIII The ground state turns out to be uncorre- 
lated (a Slater determinant) and excited states can be 
expanded in particle-hole states. We derive simple ex- 
pressions for the response functions in the spinful and 
spinless case which includes the exact solution for small 
momenta and is analytic for constant interaction. In the 
following Sec. [TV] also inter-band interactions are taken 
into account on the RPAE level. For both models the 
transition in the DSF between a curvature dominated 
limit akin to the non-interacting case and the limit of 
strong interaction close to the TL model are discussed. 
We also consider the possibility of power-law divergences 
in an intermediate regime. We close this article with 
concluding remarks in Sec. [V] 



II. RESPONSE FUNCTIONS 

We consider interacting electrons on a ring of length L. 
In units where h = 1 this is described by a Hamiltonian 
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where = k 2 /2m denotes the non-relativistic quadratic 
dispersion and v^m-MM the matrix elements of the two- 
body interaction, containing a factor 5k 1 +k 2 .k a +k4. ensur- 
ing momentum conversation. A possible spin index er is 
included in fc. 

The quantities of interest are, firstly, the (generally 
spin-dependent) density- density response function 10 

X°°' (x ~ x', t) ee ~i6(t)([ Pr7 (x, t), p a . (x', 0)]), (2) 

where p a {x) is the density operator and u = ±, which 
leads to a charge and spin response function in the spinful 
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model denoted by 



M = 



£ X ™' and x ^=Y J ^' X aa \ 



(3) 



respectively. Secondly, the dynamic structure fac- 
tor S(q, uj) being related to the response function by 
— 7r _1 Imx(9, uj) = S(q,uj) at T = and by the fluc- 
tuation dissipation theorem in general^ Switching to 
Fourier space and using a spectral decomposition we in- 
troduce spin-dependent particle-hole (ph) states 



z ) <7) = ^|o) = ^4 +9jffCfe)(T |o) ) 



(4) 



where |0) denotes the ground state to H, and get for the 
response function at T = 



X °°'{q^)=Y; 



(q,a\n)(n\q,a') (q,a'\n)(n\q,a) 



- [E n - E ] 



[E n — Eq 



(5) 

where uj = uj + iO and \n) denotes a complete set of 
eigenstates of H. Introducing states \s) = |+) + |— ) and 
| a) = |+) — |— ) the charge and spin response functions 
can be written as 

X (c) (q,u) = (slx^'ls) and X is) (q,uj) = (a\ X °°'\a). (6) 

These quantities shall subsequently be calculated in dif- 
ferent approximations. 



III. PURE INTRA-BAND INTERACTIONS 

A. (74-model and Tamm-Dancoff approximation 

Following Tomonaga, we assume the range of the inter- 
action to be much larger than the inter-particle distance 
or, equivalently, the Fourier components of the interac- 
tion are non-zero only for small momenta (smaller than 
some cut-off k c ), v(k) ^ for k < k c <C kp. The in- 
teraction leads to excitations around the Fermi points 
and one distinguishes two types of electrons: right (left)- 
movers corresponding to electrons with momenta close 
to +(—)kp. We introduce g-ology notation^ and de- 
note interactions involving only one sort of electrons as 
(?4-processes (intra-band interactions) and those between 
different types of electrons as 172— processes (inter-band 
interactions). Using this distinction the Hamiltonian is 
a sum of the kinetic energy Hq and interactions H2 and 
H4. Interestingly, for pure intra-band interactions the 
ground state of the interacting system is still the Fermi 
sea |FS), as can be seen by direct application of H 4 to 
|FS). 

We will now focus on pure (^-interaction and drop this 
limitation in Sec. IIV1 Within this frame we have com- 
plete knowledge of the ground state and the exact excited 
states can systematically be expanded in particle-hole 



states. Such a truncated expansion up to l-ph contribu- 
tions is widely used in nuclear physics and called Tamm- 
Dancoff approximation (TDA).— We write the ground 
state as a Hartree-Fock (HF) state |0) = |HF), which co- 
incides with |FS) for the (74-model, and consider excited 
states containing 1- and 2-p/i-contributions 



\n) = |HF) 



OLfX 



(7) 



where a, j3, . . . denote unoccupied (|fc| > kp) and fi, v, . . . 
occupied (\k\ < kp) one-particle states. For not too 
strong interactions the l-ph excitations are most impor- 
tant and higher orders yield only small corrections. An- 
other important point for estimating the approximation 
is the following. Using periodic boundary conditions the 
momentum is quantised, q — m q ■ L/(2n). Thus, e.g. for 
2-ph excitations, only those combinations of momenta for 
which a + /3~fi — v = q are involved. The overall num- 
ber of possible n-ph excitations is the number of possible 
partitions of a given number m q into a sum over n in- 
teger particle momenta, which is usually referred to as 
partition function in the mathematical literature.— For 
small m q the single and double ph excitations span the 
most important part of the Hilbert space. In particu- 



lar for 



8 only 1- and 2-p/i-excitations are possible 



and our approximation Eq. ([7]) yields the exact result 
for pure intra-band interactions. For m q — 20, which is 
the value for most of the numerical results in this paper, 
these excitations make up half of the states spanning the 
full Hilbert space and higher excitations give only small 
corrections. 



B. Results 

In a first step we consider single ph excitations 
and extend this approximation in the next para- 
graph. To make use of Eq. J5]) we have to calcu- 
late p/i-matrix elements of H — Eq, where Eq — E^ F 
since the ground state is a Slater determinant as ex- 
plained above. Using (HF| c^,c a ,(il - JS^c* c^HF) = 
(HF|4 / c Q ,[ J ff, C t CM ]|HF) ) one finds 



'H'\{H - E™)\an) = (e» F - e™)S aa ,6, 



where we introduced |a/x) = c^c^ |HF) as a short hand 

notation and HF energies ef F = e k +J2k> ^k,k';k,k> /( £ fc')- 
At T = the Fermi function /(e) is a simple step function 
Q(kF — \k\). Writing momenta and spin indices explicitly, 
this reads 

(RF\cl, a/ c kW (H ~ E* F )ci +q a c k jEF) = 

" £fe F )<W<W + \[9M - S^g 4 (k' - k)] (8) 

where the Fourier components 54 (k) of the interaction 
are assumed to be non- vanishing only for small momenta 
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\k\ < k c <C kp .The basic relation Eq. ((U can now be 
written as an equation for the matrix elements of the 
response function 



^|[ w -(e]^_ e HF )]&(U 
la 

-[.94(g) - 6**9i(k - l)])xw = <W*a 
or equivalently in form of a matrix equation 



(9) 
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)4fe' + [54('7)-.g4(fc-fc")] and 



(Mfc.fc' = 34(g)- The charge and spin response function 
are therefore given by 



<7<7' 



X 



2X± 



1 -9i(q)x± 
1 + . 94(g) X±' 



(11) 



leading to Xtda 



Xg/+ + X g ° 4 ,- with the HF response 
functions x±(g,w) = T^ tt ,[(w ± a)' 1 ]^ . For a con- 
stant interaction (74(g) = <?4(O)0(fc c — \q\) the latter are 
given by 



r \ 11 

X±{q,u) = - — - 

2q 2nvp 



In 



1 + a 4 — q =F Co 



(12) 



but in general have to be determined by diagonalisation 
of a. The dimensionless quantities 



a, 



■9*(0) | 
2ttvf ' 
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2fc F ' 



t'Fg 



will be used throughout this paper. Figure Q] shows 
the charge response function in 1-p/i-TDA for fixed 
momentum q/(2kp) = 0.1 and increasing interaction 
strength. A transition from the box-like shape in the 
non-interacting case to a single delta-peak as known from 
the TL model is apparent. The HF-continuum moves to 
higher energies with growing interaction and looses spec- 
tral weight to the plasmon peak. For constant interac- 
tion it is non- vanishing for frequencies 1+0:4 — g<w< 
1 + a 4 + q. 

For the spinless case the block matrix equation (flTJf col- 
lapses to a matrix equation and the response function is 
given by 



Xtda(9- w ) = X+{Q,u) +x_(g,u;), 



(13) 



with x±(g,w) as defined above, i.e. simply the HF re- 
sponse function if only 1-ph contributions are included. 

For the 2-p/i-TDA the procedure is conceptionally 
the same, but more complicated matrix elements like 
(a 1 '/3' 'n'v'\(H — E^ ¥ )\af3[iv) have to be calculated. Using 

(a'/3'fj,'u'\(H-E^)\a0fiu) = {a' n'v'WH^c^clc^Wm) 
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FIG. 1: Charge response function in TDA for Q4 = 0.0 
(dashed), Q4 = 0.2 (dashed-dotted) and Q4 = 0.4 (full line) 
for quadratic (74(9) = <?4(0)[1 — q 2 /k 2 ] interaction and constant 
momentum q = q/(2kp) = 0.1. The HF-continuum moves to 
higher energies with growing interaction strength and looses 
spectral weight compared to the emerging peak (inset). Cal- 
culations for m q — 2-irq/L — 50. The finite width of the peaks 
results from a finite imaginary part r\ = 2.5 • 10~ 3 which leads 
to a sum of Lorentzians instead of delta peaks as would be 
the case for 77 = 0. 
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FIG. 2: Response function for the spinless model in TDA 
including only 1-ph (circles)- and l-ph+2-ph (squares) contri- 
butions for 0:4 = 0.2 and quadratic interaction for constant 
momentum q/(2kp) = 0.1. The additional 2-pft-terms result 
in slight deviations only. Calculations for m, = 2nq/L = 20. 



the commutator and the resulting matrix elements have 
to be calculated. Some details of this computation are 
presented in the appendix, we only give numerical results 
and some general remarks in this section. The charge and 
spin response function keep the form given in Eq. l[TTj) but 
the functions x± differ, since the dimension of a is much 
larger now. Figure [2] shows a comparison of the response 
function for the spinless model with 1-ph and 1- and 2- 
ph excitations included for a 4 = 0.2 and q = 0.1. For 
small 0:4 the inclusion of 2-ph contributions only slightly 
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changes the resulting response function. Most of these 
additional excitations have very small weight and higher 
order excitations do contribute even less. 



C. Discussion 

The (74-model is quite remarkable considering the fact 
that the interacting ground state is still a HF state. Con- 
centrating on low-lying excited states we can use the 1-ph 
TDA to calculate response functions and obtain numer- 
ical exact results for small momenta. The results for 
both charge and spin response function are analytic for 
an interaction constant in momentum space, and for fi- 
nite interaction of arbitrary form only a numerically un- 
demanding matrix diagonalisation is necessary. Since we 
are interested in changes due to finite curvature we focus 
our discussion on the DSF where these effects are more 
obvious than in the real part of the response functions. 
The single-p/i excitations play the dominant role for weak 
interaction in a finite systems. Higher order excitations 
do of course exist, but their influence leads only to small 
corrections as can be seen in Figure [21 Endorsed with 
this control over the quality of our approximations we 
focus for the following discussion on these most domi- 
nant contributions only. 

For fixed momentum q one finds a transition from the 
box-like shape of the DSF in the non-interacting case to 
the dominance of a delta peak corresponding to the plas- 
mon of the TL model. As already pointed out by Pirooz- 
nia and Kopieta^ there are two relevant limits: 1) van- 
ishing curvature q — qj (2mvp) <C 1 or strong interaction 
a ^> 1 and 2) strong curvature g > 1 or vanishing in- 
teraction a <C 1. For the shape of the DSF only the 
ratio x = a/q of the two dimensionless quantities is im- 
portant, leading to a box- like shape for x <C 1 and to 
a dominant plasmon peak for x 3> 1. In the parame- 
ter regime between these simple limits a quite peculiar 
shape appears for intermediate ratios x « 1, the details 
depending on the type of interaction under considera- 
tion (for a comparison of different interaction forms for 
otherwise constant parameters cf. Figure [5] below) . The 
DSF exhibits a divergence at the upper (lower) bound- 
ary of the HF continuum for the charge (spin) response 
function. This shows up as an emerging peak growing 
in weight compared to the continuum with increasing in- 
teraction strength. Figure [3] shows a log-log plot of the 
charge and spin response function in 1-ph TDA for x = 1 
and x = 1/2, respectively. The two functions exhibit a 
power-law divergence at the respective boundary where 
the frequency is determined by the zero of the denomina- 
tor in Eq. (fTTj) . This picture breaks down for x ^> 1 and 
in particular for too strong interaction. The peak then 
separates from the HF continuum and power-law char- 
acteristics is clearly absent. Note that the necessary use 
of a finite imaginary part turns the delta function peaks 
into Lorentzians. These have a width of the order of 77, 
therefore a log-log plot is only reasonable down to that 
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FIG. 3: Charge (dashed) and spin (full line) response func- 
tion in 1-ph TDA close to the upper or lower boundary of 
the HF continuum in log-log plot for q = 0.1 and quadratic 
interaction with Q4 = 0.1 (charge response) and a?4 = 0.05 
(spin response), respectively. A power-law behaviour close to 
the respective boundary is discernible. Insets: Charge and 
spin response function. Calculations for m q = 2-wq/L — 20 
and 7/ = 5 • 10~ 3 . 



scale, since the explicit form of the Lorentzians appears 
below it. 

A similar behaviour holds in the spinless case. A tran- 
sition from the non-interacting case to a plasmon peak 
of the TL model at a frequency approaching uj^ l with 
increasing interaction takes place. For intermediate ra- 
tios of Q and CV4 a clear power-law behaviour shows up 
(cf. Figure [6] below). Pustilnik et al£- predicted a power- 
law divergence for the spinless model close to the lower 
boundary w_ of the non-interacting continuum for the 
DSF given by 



S{q,u) 



2q 



for < uj - uj- < q 2 /m (14) 



with an exponent /1 = ^[54(0) — gi(q)} oc x — a^/q for 
a short-ranged quadratic interaction 34(9) = 34(0) [1 — 
q 2 /k 2 ] and independent of the inter-band interaction <?2- 
The 'edge' of possible excitations is given by the lower 
bound of the HF continuum w5 F in our approxima- 
tion, and can be given explicitly for constant interaction 
.<?< = 5i(0)@(<7 — \k c \) where it reads vpq(l + a± — q). 
In contrast to Pereira et al^ the exponents we find are 
not consistent with the above mentioned predictions but 
are usually smaller. This is depicted in Figure |4] where 
the exponent extracted from log-log plots of the DSF 
in 1-ph TDA is shown for different values of the ratio 
x = cti/q along with the expected values according to 
Eq. lfl4|) . For very small x, i.e. closer resemblance to the 
non-interacting case, the agreement is good, but clearly 
deteriorates for x > 1. Note that this is not compensated 
by 1-ph contributions since they only tend to further de- 
crease the exponent. In this respect our results confirm 
the existence of a power-law, but with different exponents 
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FIG. 4: Exponent for the divergence of the DSF for the spin- 
less model as predicted by Pustilnik et a/.— (full line) according 
to Eq. Ill4[| and extracted from log-log plots within 1-ph TDA 
(circles). 

and only within a limited parameter regime where x is of 
order unity. 



IV. INTER- AND INTRA-BAND 
INTERACTIONS: THE FULL MODEL 



The previous section treated intra-band interactions 
only, i.e. we considered the pure (74-model. We now drop 
this limitation and take also inter-band interactions g 2 
into account. We treat this full model within linearised 
TDHF usually referred to as random phase approxima- 
tion with exchange (RPAE). This amounts to a clear 
improvement over single ph TDA. Correlations in the 
ground state are now included and the excited states con- 
tain more correlations in a form that particle-hole states 
can not only be created but also annihilated-^, or as Fet- 
ter and Walecka put it: "TDA has one and only one 
particle-hole pair present at any instant of time, whereas 
the RPA permits any number of particle-hole pairs to be 
present simultaneously. " 15 We first introduce the formal- 
ism, then derive our results for the charge and spin re- 
sponse function and close this section with a discussion 
of the latter. 



A. Formalism 



The change in particle-hole expectation values 5(c 



k a ,er k a +q,cr i 



in linear order due to an external time-dependent 



potential having Fourier components V£ , k (uj) — V^(q,uj)/L can be written^ 



H C k a .a C k a +q,cr) — J 



1 m*) - m: +q ) 



L U ~ y e k„+a ~ e k„ ) 



k" a" 



[54(9) - <W'540a - K)}S{cl,, a „C k ,, +q:a „)y (15) 



k" cr" 



from which we get an equation for the response function matrix elements Xkk' (<?> w ) 
matrix equation using the projection technique it reads 



uj — a 


-bi 


-b 2 


-b 2 


-64 


uj — a 


-b 2 


-b 2 


-b 2 


-b 2 


—uj — a 


-64 


-hi 


-62 


-bi 


— uj — a 



f Xpp Xpp Xpc^ Xpq \ 

Xpp Xpp Xpq Xpq 

Xqp Xqp Xqq Xqq 

V Xqp Xqp Xqq Xqq ) 



_ d5 ^ C l,^k + q ,a) 17 

_ 9v k''+ q ,k' "~ 



diag(l). 



In form of a 



(16) 



Introducing matrices (A) k a k , — (e™ q — f i k F )5kk'boo' — in Sec. [HI We thus have 

[54(g) - <W.g 4 (fc - k")] and {B 2 )% a k , = g 2 {q) the last 
equation takes a simpler form 



uj - A -B 2 
-B 2 -uj- A 



Xpp Xpq 
Xqp XQQ 



1 
1 



(17) 



Xif = (s\Xij\s) and xfj' = ( a \Xij\a) (18) 



To calculate the charge and spin response function we where i,j equal P or Q. The block matrix B 2 can ac- 
use the states |s) = |+) + |— ) and \a) = |+) — |— ) defined cordingly be written as B 2 — g 2 (q)\s)(s\. 
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B. Results 

For the charge response function Eq. fTTj) can be solved 
easily. We get two systems of two coupled equations read- 
ing e.g. 



(to - A)' 



[1 + (w - A)- 1 5 2 (o; + A)- 1 B 2 ] xpp. 



Multiplication from left with (s\ and from right with \s) 
yields with the simple structure of B 2 the PP-part (and 
QP-part as well) of the RPAE charge response. Noting 
that 



(8\(U±AY 



2Xd 



i - 94.{q)x± 

the RPAE charge response reads 



(c) 

V 



+ x£!+ + 2ft,(ff)xS + x£ ) ,- 



1 - MqWxit 



v 

A 34,- 



(19) 



(20) 



These expressions are analytic for a constant interaction 
and in general only the functions x±(q? UJ ) have to be de- 
termined numerically. Note that g 2 - and (74-interaction 
enter expression (|20"j) independently and can therefore be 
treated completely independently. Figure [5] shows the 
charge response function in RPAE for the galilei-invariant 
model where 32(9) = gi(q) = v(q) for various types of in- 
teraction, specifically for constant v q = voQ(q — \k c \), 
quadratic v q — Vo(l ~q 2 /k 2 ) and exponential interaction 
Vq = Vo exp(— \q/k c \). Note that these different forms 
of the interaction in momentum space do not lead to 
qualitatively different results, but the generic shape and 
behaviour of the DSF is similar for all of them. For pure 
intra-band interaction, i.e. g 2 = we recover the result 
of Sec. EH 

For spinless electrons the matrix equation (fT6|) is consid- 
erably simplified by dropping the spin indices and one 
finds for the RPAE response function* 



Xrpae(<?,w) 



X- +X+ + 2g 2 (g)x+X- 
1 - [92{q)] 2 X+X- 



(21) 



where x± (<?,£*;) = tEu-[( w ± a ) %k> ■ 

For the spin response function the components of the 
RPAE spin response function are needed, i.e. the func- 
tions (a\xij\a) where i,j equal P or Q. Due to the 

(s) 

symmetry of A the non-diagonal elements of Xrpae van " 
ish as the following reasoning shows. Consider \pq — 
(iomega — A) B 2 XQQ and multiply by (o| from left and 
I a) from right (analogous to the charge response), then 

HxpqH = -92(q)(a\(ui - A)" 1 |s)(s|xQQ|a) 

holds. The coefficient on the rhs vanishes, since the two 
diagonal elements to — a of to — A are equal and A is sym- 
metric. The same holds for the inverse and the expecta- 
tion value (a I ■ \s) is nothing but the difference between 
the diagonal and off-diagonal elements ^ 




1,2 1,3 

CO / (v p q) 

FIG. 5: Charge response function in RPAE for constant 



line) 



= v <d(q — |fc c |), quadratic (dashed) 



(full 

«o(l — q 2 /kZ) and exponential (dashed-dotted) interaction 
v q — vo exp( — \q/k c \). Dimensionless interaction strength 
a — 0.25 and momentum q = q/(2k,F) = 0.1 which corre- 
sponds to an intermediate relation between interaction and 
curvature of x = 2.5. Calculations for m q — 50 and fi- 
nite width of the peaks results from a finite imaginary part 
77 = 2.5 • 10" 3 . 



The RPAE spin response is thus block diagonal and 
only 



u — A 

-B 2 



-B 2 
-to — A 



\ f Xpp 


> 


- i 




H 


XQQ J 







has to be considered, yielding the simple solution 

2X- 2 X+ 



Xrpae (9,^) = J 



+ 9i(q)x- i + S4(?)x+ : 



(22) 



(23) 



which coincides with the 1-ph TDA spin response func- 
tion. This result is again analytic for a constant interac- 
tion (and coincides with earlier results- in this case) and 
easy to calculate numerically for other types of interac- 
tion. Note that the spin response function is completely 
independent of inter-band interactions, i.e. independent 
of g 2 and therefore has to be the same as in Sec. IIIII For 
increasing ratio of the interaction strength a over the di- 
mensionless momentum q a collective spin mode emerges 
which carries most of the spectral weight and moves to 
the frequency predicted by the TL model toJ)j = VFq- 



C. Discussion 

The full model shows a similar behaviour as the 34- 
model: In varying the ratio x = a/q of dimensionless 
momentum q and interaction strength a 2 u\ we find a 
transition from a box-like shape (for x S> 1 akin to the 
non-interacting case) and a dominant plasmon peak (for 
x <C 1 close to the predictions of the TL model). For 
growing interaction strength the peak approaches the fre- 
quency given by the TL model: [wjc] 2 = (vf<?) 2 [(1 + 
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2a 4 ) 2 — (2a 2 ) 2 ] for the charge mode and [wjs] 2 = (vf<?) 2 
for the spin mode. Again we find strong indications of 
power-law behaviour for the charge and spin response 
function for intermediate ratios x « 1, now not only de- 
pending on the form of the interaction under considera- 
tion but also on the strength of the inter-band interaction 
<72(0). This holds for the spinless model, too. Figured 
shows the DSF for the spinless 54-model and the full spin- 
less model for constant momentum q = 0.1. The inclu- 
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FIG. 6: Dynamic structure factor for the spinless model in 
1-ph-TBA (dashed) for a 2 = 0.1, 04 = 0.0 and RPAE (full 
line) for 02 = oa = 0.1, both for q/(2kp) = 0.1. Calculations 
for m q = 100 and 77 = 5 ■ 10" 3 . 



the divergence goes over to the dominant plasmon peak 
for values of x smaller than one. 



V. SUMMARY 

We have discussed charge and spin response functions 
for interacting electrons in one spatial dimension. We 
went beyond the assumption of the TL model and in- 
cluded deviations from a strictly linear energy disper- 
sion. We first focused on pure intra-band interactions 
(54-processes) , where the ground state is known to be a 
simple Slater determinant and expanded the excited state 
in ph excitations. Expressions for the response functions 
were derived and the DSF has been discussed. A clear 
transition from the non-interacting case to the TL model 
was found depending on the ratio of the dimensionless 
momentum and the interaction strength. For interme- 
diate values indications of power-law divergences were 
found but with exponents deviating from existing pre- 
dictions. 

The inclusion of inter-band interaction in the full model 
was realised within the RPAE and leads to quantitative 
but no qualitative changes. Again our examination of 
power-law exponents contradicted existing suggestions in 
the literature. 
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APPENDIX: TDA MATRIX ELEMENTS 



The matrix elements for the 2-ph TDA are 



(24) 



The commutators are easily calculated: [T, c^ctc^d,] = (e^ F + e§ F — e*p" — e^ F )cJ,ctc M c l/ and 



[V, cj.cjtyc] = -{ u 

mnm' 



hi a.m' a^jn^n^rrc/ ^ ^ ^mfi,m' n' ^a^ln^n' 



^mn,m' f3^rn^n^ rn ' ^ u 



Et t 1 



(25) 
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We choose a > (5 and (X > v (and similarly for other momenta) and find after some lengthy calculation the matrix 
elements 

{a' f3' p!v'\{H - Ef F )\a(3pv) = (e™ + ejp - e^ F - e^ F )5 aa >8pp'5 tltl '8y V i + v a 'p' <a p5^'5 U u' + ^ v 5 aa i5pp: 

— Sua' lyp'frPn'Svv' — Vp>n,pu>Svn' + Vp'u,Pu'^iin' ~ V p' v,Pn' & nv'~\ 
+5 a p' [Vc'^pn'Svv' — V a '^p v 'S VIM > + Va'^pv'Sfj,^' — Vi/^/i'^'] 
3pP' \ya' fj,,afj,' &uv' V a ' ^ au i 5 U ~\~ Va' u,au' ^fifi' ^a'u,afj,'^iJ,v , '\ 
+8p a i [vp'^a/M'Sw' — Vp'^av'dvfj,' + Vp'^au'S^' — W^'r/.a/i'^i/'] • 

Moreover some mixed terms of the form {\ph\(H — E^ F )\2ph) and (2ph\(H — E^ F )\lph) have to be calculated. Since 
the commutators are known already this can be done easily and one obtains 

(a'P'p,'v'\(H - E^ F )\ap) = u^g/^y/iW - v m ',n'u'S a p' + Vp' a ', av 'S^ - vp' a i iCt ^S^ 

(a' p!\{H - Ef F )\aPfXf) = V^'p^Saa' ~ Vp> a ^ v 5 a 'p + Vp a ,a>>v$w' ~ Vpa.a'^vpi'- 

With these matrix elements the response functions can be calculated numerically including 1-ph- and 2-ph- 
contributions. , 
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